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Abstract 

We have calculated the full one-loop electroweak (EW) and QCD corrections to the 
third generation scalar-fermion pair production processes e^e“ ^ 77 ^ fifiif = L t) 
at an electron-positron linear collider(LC) in the minimal supersymmetric standard model 
(MSSM). We analyze the dependence of the radiative corrections on the parameters such 
as the colliding energy ^/l and the SUSY fundamental parameters Af, tan/3, fj,, Msusy 
and so forth. The numerical results show that the EW corrections to the squark-, stau-pair 
production processes and QCD corrections to the squark-pair production processes give 
substantial contributions in some parameter space. The EW relative corrections to squark- 
pair production processes can be comparable with QCD corrections at high energies. 
Therefore, these EW and QCD corrections cannot be neglected in precise measurement 
of sfermion pair productions via 77 collision at future linear colliders. 
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1 Introduction 


The standard model (SM) has been successful in describing the strong, weak and electromag¬ 
netic interaction phenomena at the energy scale up to 10^ GeV. At the higher energy scale, 
it is likely that the minimal supersymmetric standard model (MSSM) is the most attractive 
candidate among various extensions of the SM. In the MSSM, the existence of scalar partners 
of all fermions in the SM, namely, two chiral scalar fermions /l and /r are required. At 
future colliders running in TeV energy region, the supersymmetric scalar particle // pair 
production processes are very promising channels in probing directly the existence of these 
scalar fermions, since their production cross sections can be comparatively large, if the scalar 
fermions are not too heavy. 

The two chiral SUSY states /l and fa of each fermion / turn to their mass eigenstates 
by mixing with each other. The mixing size is proportional to the mass of the corresponding 
SM fermion Generally, people believe that the sfermions of the third generation are more 
important in direct SUSY discovery than those of the former two generations, because the 
sfermions fi and fji of the third generation mix strongly to form the two mass eigenstates 
fi and f 2 - We assume that the mass eigenstates /i(/ = t,b,T) have lower masses than / 2 . 
Therefore, /i is very probably to be discovered in a relatively lower colliding energy range. 
Another significance of the sfermion pair production is that it gives access to one of the SUSY 
fundamental parameters Af, the trilinear coupling parameter. 

The future higher energy e'’'e“ linear colliders (LC) is designed to look for the evidences of 
Higgs boson and other new particles beyond the SM. There have been already some detailed 
designs of linear colliders, such as NLC m, JLC[S|, TESLAg and GLIGg. Because of the 
cleaner background of e'’'e“ collision than pp{p) collision, LG can produce more distinctive 
experimental signature of new physics. The slepton pair production at LG are intensively 
discussed in Refs. 0 13 El El- The squark pair produced by e'’'e annihilation has been 
studied thoroughly, both at tree level and at next-to-leading order ^ 3 - m the 

QCD correction to stop pair production via 77 fusion at e"''e“ linear collider is investigated. 
The scalar fermion pair production via e"’“e“ collisions e''"e“ ^ fifjif = 'T,t,b, i,j = 1,2) 
at one-loop level, has been studied in detail in usunj. They have considered the complete 
SUSY-QCD and electroweak (EW) one-loop corrections. Their results show that at the energy 
of ^/s = 500 ~ 1000 GeV, the QGD corrections are dominated while the EW corrections are 
of the same magnitude as the SUSY-QCD corrections at the higher energy scale. 

However, the future e^e“ linear colliders are designed to give other facilities operating in 
e~ e ~, 77 and other collision modes at the energy of 500 ~ 5000 GeV with a luminosity of the 
order 10 ^^ cm' -VI m The future LG’s can turn the high energy electron-positron beams 
into the Compton backscattering energetic photon beams with high efficiency in the scattering 
of intense laser photons. With the help of the new experimental techniques, it is feasible to 
yield a scaler fermion pair production directly via the high energy photon collision. Different 
options of the colliding mode are complementary to each other and will add essential new 
information to that obtained from the CERN Large Hadron Collider (LHC). Therefore, the 
sfermion pair production via 77 fusion provides another important mechanism in producing 
sfermion pair. Moreover, their production rates should be larger than those by the e'’'e“ 
annihilation because of the existing of the s-channel suppression in the latter. At the tree- 
level, the two final sfermions produced in 77 collisions should be the same sfermion mass 


2 


eigenstate, since only the electromagnetic interaction is involved. Although there are some 
studies on the e+e“ —> 77 ^ fifiif = t,b,T, i = 1, 2) at tree level|16j. the complete one-loop 
level effects of the EW and QCD in the sfermion pair production via 77 collisions are still 
absent at present. In a word, the process of scalar fermion pair production via photon-photon 
collisions e"''e“ ^ 77 ^ f^f- (/ = t,b,T,i = 1,2) will be worthwhile to investigate precisely 
and can be accessible in accurate experiments. 

In this paper, we will calculate the full one-loop EW and QCD corrections to this process. 
The paper is organized as follows: In Section 2, we give the definitions of the notations and 
the analytical calculations of the cross sections involving the 0{aew) EW and 0{as) QCD 
corrections. The numerical results and discussions are presented in Section 3. Einally, we 
give a short summery in Section 4. 


2 Analytical calculations 

In this section, we present the analytical calculations for the subprocesses 77 ^ fifiif — 
T,t,b,i = 1,2) and their parent processes e"''e“ —> 77 —> /j/j at the lowest order and the 
one-loop level in the MSSM. We adopt the ’t Hooft-Eeynman gauge and the definitions of 
one-loop integral functions in Ref. m- As we know that for the subprocesses 77 —> qiqi{q = 
t,b,i = 1, 2) there exist both QCD and EW quantum corrections, while for 77 ^ = 1, 2) 

subprocesses they have only EW quantum contributions. 


2.1 The sfermion sector and the lowest order cross section of the subpro¬ 
cess 77 ^ fifiif = i = 1, 2) 

In the MSSM, the Lagrangian mass term of the scalar fermion / can be written as 


nraass 


ifi h)^) 



if = T,t,b), 


( 2 . 1 ) 


where AIj is the mass matrix of /, expressed as 



rrifaj 

af^nuf m?- 


( 2 . 2 ) 


and 


+ifif^ — QfSiv? 0w)cos2l3m% 
rrij^ = + Qf sin^ 9w cos 2Prn% + mj, 

Qf = Ay — |u(tan/3)~^^/ . (2.3) 


where and are the soft SUSY breaking masses, is the third com¬ 

ponent of the weak isospin of the fermion, Qf the electric charge of the scalar fermion, 9\y 
the Weinberg angle, and Aj is the trilinear scalar coupling parameters of Higgs boson with 
scalar quarks, // the higgsino mass parameter. 
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The mass matrix M ? can be diagonalized by introducing an unitary matrix TZ^. The 


mass eigenstates /i, /2 are defined as 




II 

Jr 


cos 9 j sin 
— sin 6 j cos 6 j 


Then the mass term of sfermion / can be expressed 

_^rnass ^ 

where 


fi 

/2 


h 

Jr 


if = r,t,b), 




mZ 0 

ji 

0 m?~ 

J2 


The masses of fi , /2 and the angle 6j are fixed by the following equation 


(m| ,m|) = ^ [(m| - + A\af\^mj]^^^}, 


h' ' h 


2 ^ ' fh Ir fh Jr' 

2\af\mf 


tan 20 j= 2 


m^- - m j 

JL JR 


{0 < 9f < vr) 


We denote the subprocess 77 —> fifi as 

7 (Pi)+ 7 (^ 2 ) ^/i(P3) +/i(P4) if = T,t,b, i = l,2). 


(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 


where Pi andp 2 represent the four-momenta of the two incoming photons, pa andp 4 denote the 
four-momenta of the outgoing scalar fermion and its anti-particle, respectively. The momenta 
Pi{i = 1, • • • ,4) obey the on-shell equations, namely, pf = p| = 0 and p\ = p\ = . There 

Ji 

are three Feynman diagrams for this subprocess at the tree level, wiiich are shown in FiglT] 
The corresponding tree level amplitudes of this subprocess 77 —> fifi are represented as 



/z 


f^ 


h 



h/ 

/■t >< /. 


Figure 1: The lowest order diagrams for the 77 ^ fifi (/ = T,t,b) subprocess. 


Mo = M^‘ + M(f + m| 


( 2 . 10 ) 


where Mg, Mg and Mg represent the amplitudes of the t-channel, u-channel and quartic 
coupling diagrams respectively. The explicit expressions can be written as 


M* = 


- - 5 -' 

t — m~ 

Ji 


iiPl)<^RiP2) P3 pI, 


( 2 . 11 ) 
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m! = Mi (pi ^ P2), m| = 2ie^Q)g^^''e^{p^)e,{p2). (2.12) 

The Mandelstam variables t, u and s are defined as t = {pi — ps)^, it = [pi —pi)'^, s = 
{pi +P 2 )^ = {P3 +^ 4 )^. rrij, (i = 1 , 2 ) denotes the masses of the mass eigenstates of scalar 
fermions. 

The cross section at tree-level can be expressed as 


<5-0 (s) 


1 

IGvrs^ 



dt S|Mop, 


(2.13) 


with 


^max.min 


i ( 2 m| — s) ± ^/( 2 m| — s)^ — 4m^ 
2-*'^ V 


(2.14) 


The summation is taken over the spins and colors of initial and final states, and the bar over 
the summation recalls averaging over the initial spins. After integration of Eq. (ItT31) we get 
the analytical expressions of the cross section of 77 ^ /j/j subprocess at the tree level as 


^o(s) 


2710^ 

s 


Qj nI p {I + 


1-/32 


4/ x 2(1-2^2) 


(2.15) 


Here = mP- /§ and /? = — Am?- /s is the velocity of the produced scalar fermion. The 

kinematical variable v is defined as n = (1 — /3)/(l + /3). For squarks, we have = 3, while 
for sleptons = 1 . 


2.2 0 {aew) EW corrections to subprocess 77 —>■ fifi (/ = T,t,b, i = 1,2) 

In the calculation of the one-loop EW corrections, we adopt the dimensional reduction {DR) 
regularization scheme, which is supersymmetric invariant at least at one-loop level. We as¬ 
sume that there is no quark mixing, i.e., the CKM-matrix is identity matrix, and use the 
complete on-mass-shell (COMS) renormalization scheme |19j . We use FeynArts 3 m pack¬ 
age to generate the 0{aew) Feynman diagrams and amplitudes of the 0{aew) EW virtual 
contributions to 77 —> fifi (/ = r, t, b) subprocess. There are total 469 EW one-loop Feyn¬ 
man diagrams, and we classified them into four groups: self-energy, vertex, box diagrams and 
counter-term diagrams. The relevant renormalization constants are defined as 


60 — (1 + dZffje, -\- 5rnf^, 

. 1 


™z,o ~ Aq — —5ZazZ -|- (1 -|- —5Zaa)-^ 

= /i,o = 4i'7i + ^S'72, ho = zlTf2 + zll'"h. (2.16) 


where 


Z(i'^ = 6ij + -5Z(.. 

tj -^ 2 ^ 


With the on-mass-shell conditions, we can obtain the renormalized constants expressed as 

= ReY^ (m^), = ReY^^{m‘^), 
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^'ZiAA — — 


dp2 


/?eSP( 0 ) 


2 = 0 ) SZzA = 2 


mr 


= 


1 sin6lvi/i?eE|^(0) 

“2'*^^-'* + w = 2 '^* ap2 + cos%- ml ■ 

*n| =iJeSf,(m,|), Szi =-Re?A§-\, ..-, 


Ip 2 =m| ) 


24 k ) 


K = ^ 

‘j mir — m 


(z,j = l,2 ZT^i). 


/i 


( 2 . 17 ) 


( 2 . 18 ) 


Re means taking the real part of the loop integrals appearing in the self-energy. The 
= 1)2) appeared in Eos. (12.181) denote the unrenormalized sfermion self-energy in¬ 
volving only the EW interactions. The 0{aew) one-loop virtual corrections to 77 —> fifi is 
represented as 

1 pimax - 

/ di 2ReY,{M^AMl). ( 2 . 19 ) 

^min 

The expressions of tmax,min have been presented in Eq. (imi) and the summation with bar 
over head represents same operation as that appeared in Eq. (EI31)- MyA renormalized 
amplitude of the EW one-loop Eeynman diagrams, which include self-energy, vertex, box and 
counter-term diagrams. 

The 0{aew) virtual corrections contain both ultraviolet (UV) and infrared (IR) diver¬ 
gences. After renormalization procedure, the UV divergence should vanish. We have checked 
the cancellation of the UV divergence both analytically and numerically, and confirmed that 
we got a UV finite amplitude at the 0{aew) order. The IR singularity in the ig orig¬ 

inated from virtual photonic loop correction. It can be cancelled by the contribution of the 
real photon emission process. We denote the real photon emission as 


l{Pi) + 7 ( 7 * 2 ) ^ fi{pz) + fi{P 4 ) + l{k) if = T,t,b, i = 1 , 2 ), ( 2 . 20 ) 

where k = (fc®, k) is the four momentum of the radiated photon, and pi, 7 ' 2 : 7 * 3 : 7*4 are the four 
momenta of two initial photons and final state particles fifi, respectively. The real photon 
emission Feynman diagrams for the process 77 —> /i/i 7 are displayed in FigEl In our paper, 
we adopt the general phase-space-slicing method m to separate the soft photon emission 
singularity from the real photon emission process. By using this method, the bremsstrahlung 
phase space is divided into singular and non-singular regions. Then the correction of the real 
photon emission is broken down into corresponding soft and hard terms 


+ Aa^Z = + O)- (2.21) 

In the c.m.s. frame, the radiated photon energy k^ = ^\k\‘^ + is called ‘soft’ if k^ < AEj 
or ‘hard’ if k^ > AE.y. Here, m.y is a small photon mass, which is used to regulate the 
infrared divergences existing in the soft term. Although both Aaf^^ and Aa^^d depend on 
the soft photon cutoff AE^/E/j, where Ef, = ^ is the electron beam energy in the c.m.s. 
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Figure 2: The real photon emission diagrams for the process 77 ^ fifn (/ = t, t, b) 


frame, the real correction is cutoff independent. In the calculation of soft term, we use 

the soft photon approximation. Since the diagrams in FigEl with real photon radiation from 
the internal sfermion line or photon-sfermions vertex do not lead to IR-singularity, we can 
neglect them in the calculation of soft photon emission subprocesses by using the soft 

photon approximation method. In this approach the contribution of the soft photon emission 
subprocess is expressed as |19[ W]\ 


dAd 


EW 

soft 


-ddo 



27r2 


r d^k 

P3 

Pi 


_P3 • k 

Pi ■ k_ 


( 2 . 22 ) 


where the soft photon cutoff AE'.y satishes < AE^ <C \/I- The integral over the soft photon 
phase space has been implemented in Ref. d, then one can obtain the analytical result of 
the soft real photon emission correction to 77 —> fifi. 

As mentioned above, the IR divergence of the virtual photonic corrections can be exactly 
cancelled by that of soft real correction. Therefore, the sum of the 0{aew) virtual 

and soft contributions, is independent of the IR regulator m^. In the following numerical 
calculations, we have checked the cancellation of IR divergencies and verified that the total 
contributions of soft photon emission and the virtual corrections are numerically independent 
of m^. In addition, we present the numerical verification of that the total one-loop level 
EW correction to the cross section of 77 ^ /*/*, defined as Ad®^ = Aa^^ + AafJ^, is 
independent of the cutoff AE^. 
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Finally, we get an UV and IR finite O(aew) EW correction 

Act = + Aa^^^i = aod 

where 6^^ = 5^^ + <5^^ + is the 0{aew) EW relative correction. 

2.3 0{as) QCD correction to subprocess 77 —>• qiQi {q = t,b, i = 1, 2) 

In this subsection, we calculate the supersymmetric 0{as) QCD corrections. The relevant 
Feynman diagrams and the corresponding amplitudes of the subprocess 77 qiQi, {q = 
t,b, i = 1 , 2 ) both at tree-level and at one-loop level, are again generated by the package 
FeynArts 3|ini. The Feynman diagrams of the one-loop 0(as) QCD corrections also can be 
classified into self-energy, vertex, box and counter-term diagrams. The relevant renormalized 
constants used in the calculation are similar with those in the calculation of the one-loop 
0{aew) EW correction, which are defined and expressed as in Eqs. and Eqs. (ITTHl) 

respectively, except all the EW one-loop self-energies are replaced by the corresponding QCD 
ones. The SUSY QCD unrenormalized self-energy of the scalar quark g, (q = t,b,i = 1,2) 
can be written as a summation of three parts as the follows: 

4(P') = + eF(p') + (2.23) 

where denote the scalar quark self-energy parts corresponding to the 

diagrams with virtual gluon, virtual gluino exchanges and the squark quartic interactions 
respectively. The squark quartic interactions are introduced by the superpotential of the 
SUSY model. The three parts from the squark qi self-energy can be written explicitly as 

^ii^\p‘^) = {^oimg] +Ap^{Bo + Bl)\p^,mg,mg^] -h m- 9 , m-gj - Ao[mq^]} , 

^ (2.24) 

2 

^ii^\p‘^) = + {ml + mqnig sm2eq)Bo\p^,mg,mq] + p^Bi\p‘^,mg,mq]}, 

(2.25) 

eJ''^(p^) = cos'^‘^0g + sin‘^2eq}, (2.26) 

where nig denotes the small gluon mass, D = 4 — 2e is the space-time dimension and the 
group Casimir operator has Cp = ^. 

The one-loop 0(as) QCD virtual correction of subprocess 77 —> qiqi {q = t,b, i = 1,2) 
can be expressed as 

AAT® = (2.27) 

imin 

where is the renormalized amplitude of the one-loop 0{as) QCD Feynman diagrams, 

which include self-energy, vertex, box and counter-term diagrams. 

The virtual QCD corrections contain both ultraviolet (UV) and infrared (IR) diver¬ 
gences in general. To regularize the UV divergences in loop integrals, we adopt the di¬ 
mensional regularization in which the dimensions of spinor and space-time manifolds are 
extended to D = 4 — 2e. We have verified the cancellation of the UV divergence both 
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analytically and numerically. Then we get a UV finite amplitude including 0{as) virtual 
radiative corrections. The IR divergence of the QCD virtual corrections of the subprocess 
77 QiQi {q = i = 1 , 2 ) coming from virtual gluonic correction can be cancelled by the 
real soft gluonic bremsstrahlung, which is analogous to the real soft photonic one. The real 
gluon emission diagrams of the process 7 + 7 —> Qiqig are shown in Fig|21 We denote the real 
gluon emission as 

7(pi) + 7 (^ 2 ) ^ Qiips) + qi{Pi) + g{k), {q = t,b, i = 1, 2). (2.28) 



Figure 3: The real gluon emission diagrams for the process 77 —> qiqig {q = t,b, i = 1,2) 


Analogously, we use again the general phase-space-slicing method to separate the soft 
gluon emission singularity from the real gluon emission process. Therefore, the correction of 
the real gluon emission is divided into soft and hard terms 


+ ^^hard = 


soft 


hard 


soft 


'^hard 


(2.29) 


By using the soft gluon approximation, we get the contribution of the soft gluon emission 
sunbprocess expressed as 


dAd 


QCD 

soft 


-ddo 


DsCp 

2tt'^ 


r d^k 

CO 

4^ 

l\k\<AEg 

P3 -k PA-k_ 


(2.30) 


in which ^Eg is the energy cutoff of the soft gluon and < ^Eg <C V§. k^ = y\k\‘^ + rri^ 
is the gluon energy, and p 4 are the four momenta of two final state particles qi and qi. 
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In this approach, we may again refer to Ref. m to get the analytical expression of the soft 
gluon correction. 

Finally we obtain an UV and IR finite 0{as) QCD correction to the subprocess 

77 QiQi containing one-loop 0(as) QCD correction 

^ ^^QCD ^ -^jeCD (2 31) 

where + ^hard 0{as) QCD relative correction. 


2.4 The cross sections of parent processes e+e ^ 77 ^ fifi (/ = T,t,b, i = 

1 , 2 ) 

The fifi pair production via photon-photon fusion is only a subprocess of the parent process 
e'’'e“ —> 77 —> fifi. The laser back-scattering technique on electron beam can transform 
e'’'e“ beams into photon beams [SHlEllEn]. After integrating over the photon luminosity in 
an e'’“e“ linear collider, we obtain the total cross section of the process e'''e“ ^ 77 ^ fifi 
expressed 

p^max (^L ~ ~~ 

cj(s) = / dz—f^ai'j'y fifi at s = z^s), (2.32) 

72m? /y7 

Ji 

where ^/s and V§ are the e'’'e“ and 77 c.m.s. energies respectively and dL^^/dz is the 
distribution function of photon luminosity, which is expressed as 


rl T r^max 

-^ = 2z _/^/^(x)/^/,(zVx) (2.33) 

Caj^ j /yO, !rf% mJO 

^ ^ j u^rnax 

where /.^/g is the photon structure function of the electron beam |18[ I27j . For the initial 
unpolarized electrons and laser photon beams, the photon structure function is given by the 
most promising Compton backscattering as [m 1111122 ] 


/ 7 /e 


1 

m 


1 

1-x + - - 

1 — X 


Ax 4x^ 

^(1 — x) ^ 2(^1 _ xY 


(2.34) 


where 


D{0 



^) In (1 -I- ^) -I- - -1- 
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e 


2y/soJo 

mY 


1 

2(1 + 0"’ 


(2.35) 


rUe and '/s/2 represent the mass and energy of the electron respectively, uq is the laser-photon 
energy and x is the fraction of the energy of the incident electron carried by the backscattered 
photon. The maximum fraction of energy carried by the backscattered photon is Xmax = 
2ijJmax//s = £,/{! + /)■ In our calculations, we choose ujq to maximize the backscattered 
photon energy without spoiling the luminosity through e^e“ pair creation. Then we have 
^ = 2(1 + /2), Xraax ^ 0.83, and D{i) « 1.8397[in]- 
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3 Numerical results 


In this section, we present some numerical results for the one loop 0{as) QCD and 0{aew) 
EW corrections to subprocesses 77 —> fifi and the parent processes e'’'e“ ^ 77 —> /j/j. 
In our numerical calculation, the SM parameters are set to be as(m^) = 0.1190, me = 
0.5110998902 MeV, = 105.658357 MeV, = 1.77699 GeV, = 66 MeV, = 66 
MeV, me = 1.2 GeV, = 150 MeV, mb = 4.3 GeV, mt = 174.3 GeV, mz = 91.1876 
GeV, mw = 80.423 GeVjS^. There we use the effective values of the light quark masses 
{mu and m^) which can reproduce the hadron contribution to the shift in the fine structure 
constant aeu>(W' 2 )E 21 - We take the fine structure constant at the Z'^-pole as input parameter, 
aew{m‘^z)~^\j4S ~ 127.918[SIj. Then from Ea. (l2.17B we get the counter-term of the electric 
charge in DR scheme expressed as ISSIIMIIISI 


= 


6(47r)2 

2 


^ + U + log 




f '==1 

2 


m 


fk - 


+4^ ( A + log^ I (A + log ^ 


k=l 


Xk , 


k=l 


m 


Ht 


-22 ( A + log^ 


m 


w 


(3.1) 


where we take Xf = mz when mf < mz and xt = mt- Qf is the electric charge of (s)fermion 
and A = 2/e — 7 + log 47 r. Nq is color factor, which equal to 1 and 3 for (s)leptons and 
(s)quarks, respectively. It is obvious that there is a little discrepancy between our electric 
charge counter-term expression( En. (13.111 1 and that in subsection 3.1 of Ref. m- 

The MSSM parameters are determined by using EormGalc package with following input 
parameters jHH] : 

(1) The input parameters for MSSM Higgs sector are the GP-odd mass M^o and tan/3 
with the constraint tan f3 > 2.5. The masses of the MSSM Higgs sector are fixed by taking 
into account the significant radiative corrections. 

(2) The input parameters for the chargino and neutralino sector are the gaugino mass 
parameters Mi, M 2 and the Higgsino-mass parameter We adopt the grand unification 
theory(GUT) relation Mi = (5/3) tan^ 0 vi /-^2 for simplification jSI3] and the gluino mass mg 
is evaluated by mg = as{Q)/aewi^nz) sin^ 9wM2. 

(3) Eor the sfermion sector, we assume Mg = Mjy = M^ = M^ = = Msusy and 

take the soft trilinear couplings for sfermions q and I being equal, i.e., Aq = Ai = Af. 

Except above SM and MSSM input parameters, we should have some other parameters 
used in our numerical calculations, for example, the QCD renormalization scale Q, the IR 
regularization parameter m^{mg) and the soft cutoff AE^^g/Ef,. In our following numerical 
calculations, we take the QCD renormalization scale Q to be 2m j,, and set AE^^g/Ej, = 10“^, 
m^^g = 10“® GeV, if there is no other statement. As we know, the final results should be 
independent on IR regulator m^^g and the cutoff AE^^g/Ef,. Eor demonstration, we present 
the dependence of the 0{as) QCD corrections to 77 —> iiti{mi^ = 148 GeV) in conditions of 
\/I = 500 GeV and Setl parameters (see below) on the soft cutoff AEg/E^ in Eigll] The full, 


11 








dashed and dotted lines correspond to ^^hard total correction . 

As shown in this figure, the full 0(as) one-loop QCD correction Aa^'^^ is independent of 
the soft cutoff AEg/Ef, as AEg/Ef, running from 10“® to 10“^, although both and 

^^hard depend on cutoff strongly. 



Figure 4: The full 0(as) QCD corrections to 77 ^ fiti as a function of the soft gluon cutoff 
AEg/Eh in conditions of V§ = 500 GeV and 5etl parameters. 

In order to show and discuss the effects of the radiative corrections to the subprocess of 
77 ^ fifi quantitatively, we choose the following four typical data sets: 

Setl: tan/3 = 6 , M^o = 250 GeV, Msusy = 200 GeV, = 800 GeV, M 2 = 200 GeV and 
Af = 400 GeV. 

Then we have ^ = (185,223) GeV, = (148,340) GeV and = (146,250) 

GeV. 

Set2: tan /3 = 20, My^o = 300 GeV, Msusy = 400 GeV, /x = 1000 GeV, M 2 = 200 GeV and 
Af = -500 GeV. 

Then we have 2 = (354,446) GeV, ^ = (304, 533) GeV and ^ = (256, 508) 
GeV. 

Set3: tan/3 = 30, Mj^o = 300 GeV, Msusy = 250 GeV, /x = —200 GeV, M 2 = 800 GeV and 
Af = 250 GeV. 

Then we have 2 = (231,275) GeV, = (215,368) GeV and = (188,307) 

GeV. 
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SetA\ tan (5 = 30, M^o = 250 GeV, Msjjsy = 200 GeV, /x = 200 GeV, M 2 = 1000 GeV and 
= 300 GeV. 

Then we have 2 = (179,228) GeV, rni^^ = (131,346) GeV and = (124,263) 

GeV. 

With the input parameters tan/3, M^^o, Msjjsy^ A*; -^2 and Af in above data sets, all the 
masses of supersymmetric particles can be obtained by using package FormCalc. 5etl(or Set2) 
is the case of gaugino-like with small(or mediate) tan /3, but hghter(or heavier) sfermions, while 
5et3 and 5eM are higgsino-like case with larger tan/3. 

The Born and the full 0{aew) EW corrected cross sections for the 77 ^ fiTi subprocess 
as the functions of c.m.s. energy of 77 collider with above four data sets are displayed in 
Fig.5(a). There (To,i’s are the Born cross sections and di^j’s are the full O(aew) EW corrected 
cross sections for the subprocess 77 —> fjfi. The subscript i goes from 1 to 4, which correspond 
to the data Setl, Set2, Set3 and SetA respectively. The 0{a(,w) EW corrected cross section 
with SetA can achieve the maximal value 0.278 pb at the energy near the threshold \/I ~ 400 
GeV. When ^/§ approaches to 1.5 TeV, the EW corrected cross section with Set2 goes down 
to 27.5 fb, but it is still much larger than that for the process of e'’'e“ —> fifi da El with the 
same input parameters. In Fig.5(b), the relative 0{aew) EW corrections with the four data 
sets are depicted. As it can be seen in this hgure, the relative corrections 5 also have their 
maximal values at the position near the threshold energies and then decrease quantitatively 
with the increment of \/I. When the c.m.s. energy ^/§ goes from the threshold value of fifi 
pair production to 2 TeV, the full 0{aew) EW corrections can enhance or reduce the Born 
cross section depending on the colliding energy. At the position of colliding energy \/§ = 
2 TeV, the relative EW correction 5 can reach —24.6%, —24.1%, —23.5% and —23.2% for 
Setl, Set2, Set3 and SetA respectively. Fig.5(c) shows the numerical results of the cross 
sections of 77 —> f 2 f 2 subprocess both at the Born level and one-loop level, as the functions 
of the colliding energy \/I. Fig.5(d) displays the relative 0{aew) EW correction for f 2 T 2 - 
pair production as a function of VS. We find that the behavior of curves in Fig.5(c), which 
correspond to the Born, the EW corrected cross sections of f 2 'r 2 production, are quite similar 
to those in Fig.5(a) for fifi production. But the values of the cross sections for f 2 f 2 production 
are always smaller due to the heavier mass of f 2 , and can reach 0.173 pb near the threshold 
energy of f 2 f 2 pair production in the case of Setl. The magnitude of EW relative correction 
is about -26.1% or -24.5% at the position of = 2 TeV, which is close to that of 77 ^ fifi 
subprocess. 

In Fig. 6 (a) and Fig. 6 (c), we depict the full 0{aew) EW and 0(as) QGD corrected cross 
sections for the subprocess 77 —> titj. Analogously, d'o,i(* = l''')4) mean the tree-level 
cross sections corresponding to the four input data sets respectively, and di^j’s are the full 
one-loop corrected cross sections. Fig. 6 (a) demonstrates that the corresponding two curves 
for Born and 0{aew) EW corrected cross sections in the same condition of the input data set, 
have the same line shape. While Fig. 6 (c) shows obviously that the 0{as) QCD corrections 
can be larger than the 0{aew) EW corrections, especially near the threshold. The EW and 
QGD relative corrections to tjtj pair production subprocess are displayed in Fig. 6 (b) and (d)^, 
respectively. From Fig. 6 (b), we can see that the 0{aew) EW relative corrections to 77 ^ tjtj 
subprocess vary from positive values to negative ones as \/I running from the threshold value 
to 2 TeV. The absolute value of the relative correction for Set3 can reach about 34.2% 
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[GeV] ^ [GeV] 


Figure 5: (a) The Born and full O(aew) EW corrected cross sections for the 77 —> fifi 

subprocess as the functions of c.m.s. energy of 77 collider \/S with four different data sets, 
respectively, (b) The full O(aew) EW relative correction to the 77 ^ fifi subprocess. The 
solid ,dashed, dotted and dash-dotted curves correspond to four different data set cases, 
respectively, (c) The Born and full O(aew) EW corrected cross sections for the 77 ^ f 2 'f 2 
subprocess as the functions of c.m.s. energy of 77 collider \/S with four different data sets, 
respectively, (d) The full O(aew) EW relative correction to the 77 ^ f 2 'f 2 subprocess. The 
solid ,dashed, dotted and dash-dotted curves correspond to four different data set cases, 
respectively. 
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Figure 6 ; (a) The Born and full 0{aew) EW corrected cross sections for the 77 ^ iiti 

subprocess as the functions of c.m.s. energy of 77 collider \/I with four different data sets, 
respectively, (b) The full 0{aew) EW relative correction to 77 ^ titi subprocess. Four 
different curves correspond to four different data sets, respectively, (c) The Born and full 
0(as) QCD corrected cross sections for the 77 ^ titi subprocess as the functions of c.m.s. 
energy of 77 collider v/I with four different data sets, respectively, (d) The full 0(as) QCD 
relative correction to 77 ^ titi subprocess. 
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Figure 7: (a) The Born and full 0{aew) EW corrected cross sections for the 77 ^ ^ 2^2 

subprocess as the functions of c.m.s. energy of 77 collider \/I with four different data sets, 
respectively, (b) The full 0{aew) EW relative correction to 77 ^ ^ 2^2 subprocess. Four 
different curves correspond to four different data sets, respectively, (c) The Born and full 
0(as) QCD corrected cross sections for the 77 ^ ^ 2^2 subprocess as the functions of c.m.s. 
energy of 77 collider v/I with four different data sets, respectively, (d) The full 0(as) QCD 
relative correction to 77 ^ t 2 t 2 subprocess. 
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at the point of \/I = 2 TeV and be over 32% near the threshold for Setl. Furthermore, with 
the same input parameters as used in CSI, for example, Set2, our calculation shows that, when 
\/I is between 1200 GeV and 2000 GeV, the EW relative correction to 77 ^ iiti subprocess 
is about —24.1 ~ —31.8%, while the EW relative correction to e'’'e“ —> titi process is about 
— 10% We notice that on the curves in Eig. 6 (b) there are some small spikes which are 
due to the resonance effects. For example, the resonance effect at the position of ^/§ ~ 525 
GeV is caused by ~ 2mfj+ for input data Setl, while the resonance effect at the position 
of \/I ~ 1066 GeV is caused by \/| ~ for input data Set2. Furthermore, when we 
observe Fig. 6 (d), which shows the relative QGD correction as a function of c.m.s. energy \/I 
for 77 ^ iiti, we can see that the values of decrease rapidly to minimal values after 

\/I just goes up from the threshold energy and then increase slowly to 23.8%, 11.4%, 19.3% 
and 27.6% at the position of V§ = 2 TeV for Setl, Set2, Set3 and Set4, respectively. 

The results for subprocess 77 —> ^ 2^2 are drawn in Fig.7(a-d). The full 0{aew) EW 
and 0(as) QGD corrected cross sections are plotted in Fig.7(a) and Fig.7(c), respectively. 
Gomparing these two figures with Fig. 6 (a) and Fig. 6 (c), we can see that the cross sections 
for the 77 —> t 2 t 2 subprocess are almost one order smaller than those for the 77 ^ titi 
subprocess quantitatively due to The EW and QGD relative corrections to 

subprocess 77 ^ ^ 2^2 are plotted in Eig.7(b) and Eig.7(d), respectively. From Fig.7(d), we 
see that the values of QGD relative corrections are rather large near the threshold, and at 
position of \/I = 2 TeV they are 10.2%, —0.7%, 12.5% and 13.2% for data Setl, Set2, Set3 
and Set4 respectively, which are less than the corresponding QGD relative corrections to the 
77 ^ 1^1 subprocess shown in Fig. 6 (d). The absolute EW relative corrections are generally 

larger than those of the absolute QGD relative corrections shown in Fig.7(d), except in the 
threshold energy regions. The values of EW relative corrections are about —45% for the two 
gaugino-like data sets Setl and Set2 and —43% for the two higgsino-like data sets Set3 and 
Set4 at the position of \/I = 2 TeV. The EW relative corrections to 77 —> QQ cross section 
are negative in the range of = 700 ~ 2000 GeV with all the four data sets and have 
the minimal values near the position of \/I = 800 GeV for Setl, Set3 and Set4, and in the 
vicinity of \/I ~ 1200 GeV for Set2. 

We also show the bibi, {i = 1,2) pair productions in Fig .8 and Fig.9. Fig. 8 (a) is plotted 
for the Born and full 0{aew) EW corrected cross sections of the 77 —> bibi subprocess as the 
functions of V§ with four data sets respectively, and Fig. 8 (c) for the Born and full 0{as) QGD 
corrected cross sections. Fig. 8 (b) and Fig. 8 (d) display the EW and QGD relative corrections, 
respectively. These figures show that the behaviors of the curves for bibi pair production 
are similar to those for titi pair production. Comparing Eig. 6 (b) with Fig. 8 (b), we notice 
that the full 0{aew) EW relative correction to 77 —> 6161 subprocess is larger than that to 
the 77 —> titi subprocess. The maximum absolute value of the former can reach 50.7% for 
Set2 at the position of \/l ~ 2000 GeV, which is even larger than the QGD correction to 
77 • Again, all of the small spikes appearing on the curves of Eig. 8 (a-b) are due to 

the resonance effects, such as the spikes at the positions of \/I ~ 525 GeV for Setl and 
\/I ~ 621 GeV for Set2 are caused by \/I ~ 2mij+, while condition of ~ 2^-12 leads to the 
spike at the position of \/j ~ 1066 GeV for Set2. In Eig. 8 (d), the solid, dashed, dotted and 
dash-dotted lines correspond to the QGD relative corrections with parameter scenarios Setl, 
Set2, Set3 and Set4, respectively. Although Fig. 8 (d) demonstrates that the the QGD relative 
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Figure 8 ; (a) The Born and full 0{aew) EW corrected cross sections for the 77 ^ bibi 

subprocess as the functions of c.m.s. energy \/I with four data sets, respectively, (b) The 
full one-loop 0{aew) EW relative correction to 77 ^ bibi subprocess, (c) The Born and full 
0(as) QCD corrected cross sections for the 77 ^ bibi subprocess as the functions of c.m.s. 
energy \/I with four data sets, respectively, (d) The full 0(as) QCD relative correction to 
77 ^ subprocess. 
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Figure 9; (a) The Born and full O(aew) EW corrected cross sections for the 77 ^ 6262 

subprocess as the functions of c.m.s. energy \/I with four data sets, respectively, (b) The full 
O(aew) EW relative corrections to 77 ^ 62^2 subprocess, (c) The Born and full 0(as) QCD 
corrected cross sections for the 77 ^ 62^2 subprocess as the functions of c.m.s. energy a/I 
with four data sets, respectively, (d) The full 0(as) QCD relative correction to 77 ^ 6262 
subprocess. 
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corrections in the region near the threshold energy of the bibi pair production are extremely 
large, these values are untrustworthy due to the non-perturbative QCD effects. The values 
of the relative 0(as) QCD corrections at the position of y/I = 2000 GeV, are 23.1%, 13.4%, 
21.5% and 28.7%, for Setl, Set2, Set3 and Seti respectively. 

The full 0{aew) EW and 0(as) QCD corrected cross sections for the 77 ^ 62^2 subprocess 
are depicted in Fig.9(a) and Fig.9(c) separately, while their corresponding relative corrections 
are plotted in Fig.9(b) and Fig.9(d) respectively. Although the line-shapes in Fig.9(a) and 
fig.9(c) are similar with the corresponding ones in Fig. 8 (a) and (c) for the bibi pair production, 
the values of the corrected cross section in Fig.9(a) and Fig.9(c) are much smaller due to 

Nevertheless, Fig.9(b) shows that when \/I is large enough, the absolute EW 

relative corrections to 77 ^ 6262 approach about 50% or beyond for all four data sets. The 
peak at the position of \/I ~ 525 GeV for Setl in Fig.9(b) comes from the resonance effect 
of \/I 2mff+. We find also from Fig.9(d) that the absolute QCD relative corrections to 
the 77 ^ 62^2 subprocess are generally comparable with the EW corrections or even smaller 
than EW ones, especially in the large colliding energy region. 

In the following discussion, we present some numerical results about the parent process 
e+e” ^ 77 ^ fifiif = T,t,b, i = 1,2). For convenience, we denote the cross sections of the 
parent process e'''e“ ^ 77 ^ /j/j containing the 0{aew) EW and 0{as) QCD corrections as 

= ao + = uo(l + 5^^), = uq + Aa^^^ = ao(l + 5^^^) 

where 5^^ and are the 0{aew) EW and 0(as) QCD relative correction respectively. In 

the following numerical calculations, we take the input parameters of higgsino-like data Set3, 
but let Msusy running from 200 Gev to 400 GeV. Fig.10(a) and Fig.10(c) show the Born and 
full 0{aew) EW corrected cross sections for the e'''e“ ^ 77 ^ fifi and e'''e“ ^ 77 ^ f 2 r 2 
process as the functions of the soft-breaking sfermion mass Msusy- In Fig. 10 (a), the solid 
and dashed curves correspond to Born and full 0{aew) EW corrected cross sections for ^/s = 
500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively . It is obvious that all the curves for the 
Born and EW corrected cross sections decrease rapidly as Msusy going up from 200 to 400 
GeV, but the damping decrement is getting smaller with the increment of the colliding energy 
■y/s. We can read from Fig. 10(a) that the values of the EW corrected cross sections decrease 
from 28.3 fb, 98.9 fb, 110 fb and 92.8 fb to 2.25 fb, 0.1 fb, 1.28 fb and 20.8 fb for ^ = 500 
GeV, 800 GeV, 1000 GeV, 2000 GeV respectively, when Msusy increases from 200 GeV to 
400 GeV. In Fig.lO(c), the curves show the Born and 0{aew) EW corrected cross sections for 
e+e“ ^ 77 ^ f 2 f -2 process with y/s = 800 GeV, 1000 GeV, 2000 GeV respectively. All the 
curves have the analogous tendency as the curves in Fig.10(a). The values of the corrected 
cross sections in Fig. 10(c) are smaller than those for e+e“ ^ 77 ^ fiTi in Fig. 10(a). The 
EW relative corrections to the e'''e“ ^ 77 ^ fifi process as the functions of Msusy are 
depicted in Fig.lO(b) for y/s = 500 GeV, 800 GeV, 1000 GeV and 2000 GeV. From this 
figure, we can see that in the range of Msusy = 200 GeV to 400 GeV, this relative correction 
can reach —5.46% at the position of Msusy = 340 GeV when we take y/s = 800 GeV. 
If we take e'’'e“ colliding energy y/s = 2 TeV, we can get —18.69% relative correction to 
e'''e“ ^ 77 ^ fifi process when Msusy = 400 GeV. Fig.lO(d) displays the EW relative 
corrections to e'’“e“ ^ 77 ^ 72^2 process. We can see from this figure that the numerical 
values of these relative corrections increase rapidly from —16.43% to —3.56% when Msusy 
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Figure 10: (a) The Born and full O(aew) EW corrected cross sections for the e'''e“ —> 77 ^ 
fi 77 process as functions of the soft-breaking sfermion mass Msusy with ^/s = 500 GeV, 
SOOGeV, 1000 GeV, 2000 GeV, respectively, (b) The full O(aew) EW relative corrections to 
the e+e“ ^ 77 ^ fi 77 process as the functions of Msusy with ^/s = 500 GeV, 800 GeV, 
1000 GeV, 2000 GeV, respectively, (c) The Born and full 0{aew) EW corrected cross sections 
for the e+e" ^ 77 ^ f 27^2 process as functions of the soft-breaking sfermion mass Msusy 
with = 500 GeV, SOOGeV, 1000 GeV, 2000 GeV, respectively, (d) The full 0{aew) EW 
relative corrections to the e'''e“ —> 77 ^ f 2 'f 2 process as the functions of Msusy with ^/s = 
500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 
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goes up from 200 GeV to 310 GeV for y/s = 800 GeV, and from —17.16% to —11.7% when 
Msusy increases from 200 GeV to 390 GeV for y/s = 1000 GeV. But for y/s = 2000 GeV 
the value of 5^^ is relative stable and keeps in the range of [-18.93%, -19.72%] as Msusy 
varying from 200 GeV to 400 GeV. 

The numerical results for the process e'’'e“ ^ 77 ^ titi are plotted in Fig.11. Fig.11(a) 
and Fig. 11(c) display the Born and full one-loop EW and QCD corrected cross sections as 
the functions of Msusy with y/s = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 
As we expect, the curves in Fig.11(a) for the cross sections in Born approximation and at 
0{aew) EW one-loop level, have some similar behaviors with those for the fifi production 
process shown in Fig.10(a). We can find from Fig.11(a) that the full 0{aew) EW corrected 
cross sections decrease from 87.3 fb and 65.3 fb to 1.21 fb and 11.8 fb when Msusy goes 
from 200 GeV to 400 GeV for y/s = 1000 GeV, 2000 GeV, respectively. In Eigll(c) the solid 
and dotted curves correspond to the Born and one-loop QCD corrected cross sections versus 
Msusy with y/s = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. Erom this figure, 
we can see that the value of the QCD corrected cross section reaches 118 fb for y/s = 1000 
GeV with our chosen parameters. In order to study the EW and QCD radiative corrections 
more clearly, we plot the EW and QCD relative corrections to the e''“e“ ^ 77 ^ titi process 
in Eig.ll(b) and Fig.ll(d). In Fig.ll(b), the resonance effect at the position of Msusy = 386 
GeV is due to the condition of ~ mt + m^o. Fig. 11(b) shows the EW relative corrections 
for y/s = 1000 GeV and 2000 GeV can reach —27.15% and —26.78% at the position of Msusy 
= 400 GeV. Erom Eig.ll(d) we find the curves for y/s = 500 GeV, 800 GeV and 1000 GeV 
go up fleetly with the increment of Msusy 1 but for the curve of y/s = 2000 GeV the relative 
correction is almost stable, varying in the range of [14.6%, 13.3%]. 

The results for e^e“ ^ 77 ^ t 2 t 2 are represented in Fig.12. Fig.12(a) shows the plot of 
the Born and full one-loop EW corrected cross sections versus Msusy- Fig. 12(b) describes 
the EW relative corrections as the functions of Msusy- The QCD corrected cross sections 
and QCD relative corrections are plotted in Fig.12(c) and Fig.12(d), respectively. In Fig.12(a) 
the EW corrected cross section of e+e“ ^ 77 ^ ^ 2^2 process decreases from 13.2 fb to 5.3 fb 
for y/s = 2000 GeV, when Msusy increases from 200 GeV to 400 GeV. At the position of 
Msusy = 266 GeV in Fig. 12(a), there is a dithering on the curve of y/s = 2000 GeV, which is 
due to resonance effect of ^ mt + m^o. The resonance effect on the EW relative correction 

curves at the position near Msusy = 266 GeV is also shown in Fig. 12(b). In Fig. 12(d), the 
QCD relative corrections for y/s = 1000 GeV are rather larger and vary in the region between 
36.4% and 76.1% with the increment of Msusy, but the QCD relative corrections for y/s = 
2000 GeV are smaller, and have the values in the range between 15.5% and 18.6%. 

We also present the results of the e+e“ ^ 77 ^ bibi process in Fig.l3. Fig.l3(a) and 
Fig.13(c) show the Born and full one-loop EW and QCD corrected cross sections, respectively. 
In Fig.13(a) and Fig.13(c) we find that all the curves of cross sections at the Born level and 
involving EW and QCD one-loop contributions, decrease with the increment of the Msusy- 
For example, when Msusy varies from 200 to 400 GeV, the two curves of the cross sections 
including full one-loop EW corrections for y/s = 1000 GeV and 2000 GeV in Fig. 13(a), goes 
down from 11.1 fb and 6.6 fb to 0.13 fb and 0.65 fb respectively. While the two curves in 
Fig.13(c), which represent the cross sections including QCD corrections for y/s = 1000 GeV 
and 2000 GeV, decrease from 12.8 fb and 8.4 fb to 0.3 fb and 1.3 fb, respectively. Fig.l3(b) 
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Figure 11; (a) The Born and full O(aew) EW corrected cross sections for the e"''e“ ^ 

77 ^ process as the functions of the soft-breaking sfermion mass Msusy with ^/s = 
500, 800, 1000, 2000 GeV, respectively, (b) The full O(aew) EW relative corrections to the 
e+e” —> 77 ^ process as the functions of Msusy with ^/s = 500 GeV, 800 GeV, 1000 
GeV, 2000 GeV, respe_ctively. (c) The Born and full 0(as) QCD corrected cross sections for 
the e + e “ ^ 77 ^ iiti process as the functions of the soft-breaking sfermion mass Msusy 
with = 500, 800, 1000, 2000 GeV, respectively, (d) The full 0(as) QGD relative corrections 
to the e ''' e “ ^ 77 ^ titi process as the functions of Msusy with ^/s = 500 GeV, 800 GeV, 
1000 GeV, 2000 GeV, respectively. 
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Figure 12: (a) The Born and full 0{aew) EW corrected cross sections for the e'’'e“ ^ 77 ^ 
^ 2^2 process as the functions of the soft-breaking sfermion mass Msusy with ^/s = 1000 GeV, 
2000 GeV, respectively, (b) The full O(aew) EW relative corrections to the e^e“ ^ 77 ^ t 2 t 2 
process as the functions of Msusy with ^/s = 1000 GeV, 2000 GeV, respectively, (c) The Born 
and full 0(as) QGD corrected cross sections for the e'''e“ ^ 77 ^ ^ 2^2 process as functions 
of the soft-breaking sfermion mass Msusy with ^/s = 1000 GeV, 2000 GeV, respectively, (d) 
The full 0(as) QGD relative corrections to the e'''e“ —> 77 ^ ^ 2^2 process as the functions 
of Msusy with ^/s = 1000 GeV, 2000 GeV, respectively. 
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Figure 13: (a) The Born and full 0{aew) EW corrected cross sections for the e+e —> 77 ^ 
bibi process as the functions of the soft-breaking sfermion mass MguSY with ^/s = 500 GeV, 
800 GeV, 1000 GeV, 2000 GeV, respectively, (b) The full 0{aew) EW relative corrections to 
the e+e” ^ 77 ^ bibi process as the functions of Msusy with ^/s = 500 GeV, 800 GeV, 
1000 GeV, 2000 GeV, respectively, (c) The Born and full 0(as) QGD corrected cross sections 
for the e'’'e“ —> 77 ^ bibi process as the functions of the soft-breaking sfermion mass Msusy 
with = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively, (d) The full 0{as) QGD 
relative corrections to the e+e" ^ 77 ^ bibi process as the functions of Msusy with ^/s = 
500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 
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Figure 14; (a) The Born and full 0{aew) EW corrected cross sections for the e"*'e ^ 

77 ^ 62^2 process as the functions of the soft-breaking sfermion mass Msusy with ^/s = 
800 GeV, 1000 GeV, 2000 GeV, respectively, (b) The full 0{aew) EW relative corrections 
to the e+e“ ^ 77 ^ 62^2 process as the functions of Msusy with ^/s = 800 GeV, 1000 
GeV, 2000 GeV, respectively, (c) The Born and full 0(as) QCD corrected cross sections for 
the e^e“ —>■ 'jj —>■ 62^2 process as the functions of the soft-breaking sfermion mass Msusy 
with ^/s = 800 GeV, 1000 GeV, 2000 GeV, respectively, (d) The full 0{as) QGD relative 
corrections to the e+e" ^ 77 ^ 6262 process as the functions of Msusy with ^/s = 800 GeV, 
1000 GeV, 2000 GeV, respectively. 
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shows the full one-loop EW relative corrections to the e+e” ^ 77 ^ bibi process, as the 
functions of MsuSY for ^/s = 500, 800, 1000, 2000 GeV respectively. We see that there 
occur the resonance effects on each curve at the position of ^ mfj+ — ~ 133 GeV 

(corresponding to MsuSY ~ 212 GeV). Fig. 13(d) displays the full one-loop QCD relative 
corrections as the functions of Msusy with ^/s = 500, 800, 1000 , 2000 GeV, respectively. 
The QCD relative corrections can be rather larger for ^/s = 500 GeV and 800 GeV, and can 
reach 59% and 69% at the positions of Msusy = 250 GeV and 350 GeV, respectively. 

Finally, we present the Born and full one-loop EW and QCD corrections to e'’'e“ ^ 
77 ^ 62^2 process in Fig.14. Comparing Fig.14(a) with Fig.13(a), we can see that the Born 
and full one-loop EW corrected cross sections for 62^2 pair production are smaller than the 
corresponding ones for bibi pair production because of But the EW relative 

corrections to e^e“ ^ 77 ^ 62^2 process are rather large and can be comparable with their 
QCD corrections as shown in Fig.14(b) and Fig.14(d). We can also find that when Msusy 
= 200 GeV, the EW relative corrections are -57.8% (for ^/s = 800 GeV), -56.5% (for ^/s = 
1000 GeV) and -52.4% (for -y/i = 2000 GeV), and when Msusy = 200 GeV the QCD relative 
corrections are 38.1% and 24.1% for = 800 GeV and 1000 GeV respectively, but for ^/s 
= 2000 GeV the QCD relative correction varies in the a small range of [15.4%, 16.9%]. 


4 Summery 

In this paper, we have calculated the full 0{aew) EW and 0(as) QCD contributions to 
the third generation scalar fermion {fi,ti,bi,i = 1 , 2 ) pair production in 77 collision at 
an e'*“e“ collider. The calculation of the radiative corrections was carried out analytically 
and numerically. The numerical results were discussed in conditions of both gaugino-like 
and higgsino-like input parameter scenarios. Our investigation shows that the full 0{aew) 
EW relative corrections to both subprocesses and parent processes are typically of the order 
10~30%, and the EW relative corrections to the squark pair production can be compara¬ 
ble with the 0(as} QCD contributions in some parameter space, especially in high 77 and 
e'’'e“ colliding energy regions. For example, the EW relative corrections can reach -50% 
and -40.76% at V§(\/s) = 2000 GeV to the subprocess 77 ^ 6161 and its parent process 
e+e” ^ 77 ^ bibi, respectively. We find the full 0{as) QCD corrections to these squark 
pair production subprocesses are also large under our input data sets, for example, the QCD 
relative correction is 28.7% for 77 ^ 6161 subprocess with y/s = 2000 GeV and SetA pa¬ 
rameters. In conclusion, our numerical results have indicated that the full 0{aew) EW cor¬ 
rections to e'''e“ ^ 77 ^ (/ = t,b,T,i = 1,2) processes and 0{as) QCD corrections 

to e+e“ —>■ —>■ (^q = t,b,i = 1 , 2 ) processes, give substantial contributions in some 

parameter space. Therefore, these radiative corrections cannot be neglected in considering 
the third generation sfermion pair productions in 77 collision mode at future linear colliders. 
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Figure Captions 

Figure 1 The lowest order Feynman diagrams for the subprocess 77 ^ fifi (/ = r, t, b). 
Figure 2 The real photon emission diagrams for the process 77 ^ fifi'J (/ = t, t, b). 
Figure 3 The real gluon emission diagrams for the process 77 ^ qiQig {q = t,b). 

Figure 4 The full 0{as) QCD corrections to 77 — > titi as a function of the soft gluon 
cutoff AEg/Efy in conditions of V§ = 500 GeV and Setl parameters. 
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Figure 5(a) The Born and full O(aew) EW corrected cross sections for the 77 ^ firi 
subprocess as the functions of c.m.s. energy of 77 collider \/S with four different data sets, 
respectively. 

Figure 5(b)The full O(aew) EW relative correction to the 77 ^ fifi subprocess. The 
solid ,dashed, dotted and dash-dotted curves correspond to four different data set cases, 
respectively. 

Figure 5(c)The Born and full O(aew) EW corrected cross sections for the 77 ^ 7^277 
subprocess as the functions of c.m.s. energy of 77 collider \/I with four different data sets, 
respectively. 

Figure 5(d) The full O(aew) EW relative correction to the 77 ^ f 2 'f 2 subprocess. The 
solid ,dashed, dotted and dash-dotted curves correspond to four different data set cases, 
respectively. 

Figure 6(a) The Born and full 0{aew) EW corrected cross sections for the 77 ^ titi 
subprocess as the functions of c.m.s. energy of 77 collider \/I with four different data sets, 
respectively. 

Figure 6(b) The full 0{aew) EW relative correction to 77 ^ titi subprocess. Four 
different curves correspond to four different data sets, respectively. 

Figure 6 (c) The Born and full 0{as) QCD corrected cross sections for the 77 ^ titi 
subprocess as the functions of c.m.s. energy of 77 collider \/I with four different data sets, 
respectively. . 

Figure 6(d) The full 0(as) QCD relative correction to 77 ^ iiti subprocess. 

Figure 7(a) The Born and full 0{aew) EW corrected cross sections for the 77 ^ t 2 Q 
subprocess as the functions of c.m.s. energy of 77 collider \/I with four different data sets, 
respectively. 

Figure 7(b) The full 0{aew) EW relative correction to 77 ^ QQ subprocess. Four 
different curves correspond to four different data sets, respectively. 

Figure 7(c) The Born and full 0{as) QCD corrected cross sections for the 77 ^ QQ 
subprocess as the functions of c.m.s. energy of 77 collider \/I with four different data sets, 
respectively. 

Figure 7(d) The full 0(as) QCD relative correction to 77 ^ QQ subprocess. 

Figure 8(a) The Born and full 0{aew) EW corrected cross sections for the 77 ^ bibi 
subprocess as the functions of c.m.s. energy \/I with four data sets, respectively. 

Figure 8(b) The full one-loop 0{aew) EW relative correction to 77 ^ bibi subprocess. 

Figure 8(c) The Born and full 0(as) QCD corrected cross sections for the 77 ^ 6161 
subprocess as the functions of c.m.s. energy \/I with four data sets, rj 3 spectively. 

Figure 8(d) The full 0(as) QCD relative correction to 77 ^ 6161 subprocess. 

Figure 9(a) The Born and full 0{aew) EW corrected cross sections for the 77 ^ 6262 
subprocess as the functions of c.m.s. energy \/I with four data sets, re_spectively. 

Figure 9(b) The full 0{aew) EW relative corrections to 77 ^ 6262 subprocess. 

Figure 9(c) The Born and full 0{as) QCD corrected cross sections for the 77 ^ 6262 
subprocess as the functions of c.m.s. energy \/I with four data sets, respectively. 

Figure 9(d) The full 0(as) QCD relative correction to 77 ^ 6262 subprocess. 

Figure 10(a) The Born and full 0{aew) EW corrected cross sections for the e'’'e“ ^ 
77 ^ fifi process as functions of the soft-breaking sfermion mass Msjjsy with y/s = 500 
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GeV, SOOGeV, 1000 GeV, 2000 GeV, respectively. 

Figure 10(b) The full 0{aew) EW relative corrections to the e+e" ^ 77 ^ fifi process 
as the functions of Msusy with ^/s = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 10(c) he Born and full ©(ugu,) EW corrected cross sections for the e'’'e“ ^ 
77 ^ f 2 'f 2 process as functions of the soft-breaking sfermion mass Msusy with ^/s = 500 
GeV, SOOGeV, 1000 GeV, 2000 GeV, respectively. 

Figure 10(d) The full 0{aew) EW relative corrections to the e+e” ^ 77 ^ f' 2 r 2 process 
as the functions of Msusy with y/s = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 11(a) The Born and full 0{aew) EW corrected cross sections for the e'’'e“ ^ 
77 ^ iiti process as the functions of the soft-breaking sfermion mass Msusy with ^/s = 500, 
800, 1000, 2000 GeV, respectively. 

Figure 11(b) The full O(aew) EW relative corrections to the e'''e“ ^ 77 ^ titi process 
as the functions of Msusy with ^/s = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 11(c) The Born and full 0(as} QGD corrected cross sections for the e"’“e“ ^ 
77 ^ titi process as the functions of the soft-breaking sfermion mass Msusy with ^/s = 500, 
800, 1000, 2000 GeV, respectively. 

Figure 11(d) The full 0{as) QGD relative corrections to the e'''e“ ^ 77 ^ titi process 
as the functions of Msusy with = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 12(a) The Born and full 0{aew) EW corrected cross sections for the e'’'e“ ^ 
77 ^ QQ process as the functions of the soft-breaking sfermion mass Msusy with ^/s = 
1000 GeV, 2000 GeV, respectively. 

Figure 12(b) The full 0{aew) EW relative corrections to the e+e" ^ 77 ^ QQ process 
as the functions of Msusy with ^/s = 1000 GeV, 2000 GeV, respectively. 

Figure 12(c) The Born and full 0(as} QGD corrected cross sections for the e'*“e“ —> 
77 ^ QQ process as the functions of the soft-breaking sfermion mass Msusy with ^/s = 
1000 GeV, 2000 GeV, respectively. 

Figure 12(d) The full 0{as) QGD relative corrections to the e'''e“ —>■ —>■ QQ process 

as the functions of Msusy with ^/s = 1000 GeV, 2000 GeV, respectively. 

Figure 13(a) The Born and full 0{aew) EW corrected cross sections for the e'’'e“ —> 
77 ^ process as the functions of the soft-breaking sfermion mass Msusy with ^/s = 500 
GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 13(b) The full 0{aew) EW relative corrections to the e"’'e“ ^ 77 ^ 6161 process 
as the functions of Msusy with -y/s = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 13(c) The Born and full 0(as} QGD corrected cross sections for the e'’“e“ ^ 
77 ^ process as the functions of the soft-breaking sfermion mass Msusy with -y/i = 500 
GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 13(d) The full 0{as) QGD relative corrections to the e'’'e“ ^ 77 ^ bibi process 
as the functions of Msusy with = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 14(a) The Born and full 0{aew) EW corrected cross sections for the e'’'e“ —> 
77 ^ 6262 process as the functions of the soft-breaking sfermion mass Msusy with ^/s = 800 
GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 14(b) The full O(aew) EW relative corrections to the e"’'e“ ^ 77 ^ 62^2 process 
as the functions of Msusy with ^/s = 800 GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 14(c) The Born and full 0(as} QGD corrected cross sections for the e'*“e“ ^ 
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77 ^ ^ 2^2 process as the functions of the soft-breaking sfermion mass M^usy with ^/s = 800 
GeV, 1000 GeV, 2000 GeV, respectively. 

Figure 14(d) The full 0{as) QGD relative corrections to the e'*'e“ —> 77 ^ 62^2 process 
as the functions of Msusy with ^/s = 800 GeV, 1000 GeV, 2000 GeV, respectively. 
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